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Caveats 

The  model  described  In  this  report  was  oriptnally  Intended  for  use  In  connection  with 
nearly  vertical  laser  beam  propagation  (ground  to  space  and  vice  versa).  For  this  reason  the 
verification  was  based  on  such  trajectories  and  therefore  this  model  has  not  been  tested  for 
horizontal  or  nearly  horizontal  propagation. 

The  geometry  of  atmospheric  turbulence  has  been  described  in  the  literature  as  being 
like  “horizontal  pancakes".  The  layers  of  turbulence  are  thin  in  the  vertical  and  very  large  In 
the  horizontal.  This  means  that  the  case  for  nearly  horizontal  beam  propogation  Is  very  dif¬ 
ferent  from  the  nearly  vertical  case.  The  Cn2  values  change  rapidly  and  wildly  In  the  vertical 
direction  and  can  be  almost  uniform  horizontally  over  the  distance  of  kilometers.  What  Is 
required  is  validation  for  horizontal  trajectories,  and  the  latter  has  yet  to  be  done. 

There  Is  yet  another  point  to  be  made.  The  piesent  model  presumes  that  the  beam  is 
nearly  perpendicular  to  the  "tubulent  pancake  layers”.  What  happens  when  the  beam  is  nearly 
parallel?  It  Is  possible  that  it  could,  for  example,  reflect  from  sharp  discontinuities  of  Index 
of  refraction  at  layer  boundaries.  This  and  other  such  effects  are  not  Included  in  the  physics 
of  the  piesent  mode!.  This  may  have  to  be  Investigated  also  before  one  has  a  vaiid  approach 
to  the  case  of  nearly  horizontal  beams.  At  the  very  least  some  theoretical  justification  would 
have  to  be  given  for  Ignoring  these  omissions. 
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A  Model  for  C*  (Optical  Turbulence)  Profiles 

Using  Radiosonde  Data 

1.  Introduction 

Currently,  models  are  needed  to  convert  standard  meteorological  data  Into  veitlcal  pro¬ 
files  of  Cn2,  the  structure  constant  for  optical  Index  of  refraction  fluctuations.  For  background 
information  see  Tatarskl.1  (1961),  and  Dewan2  (1980).  The  latter  report  contains  many  other 
references.  Profiles  of  C,,2  from  the  ground  to  20  or  30  km  are  needed  to  ascertain  the  effects 
of  turbulence  on  laser  beam  propagation  from  ground  to  space  as  well  as  on  light  propaga¬ 
tion  from  space  to  ground.  Such  Information  would  be  available  In  large  quantities  tf  radio¬ 
sonde  information  could  be  converted  to  C,2  profiles.  Such  Information  could  then  be  used 
for  assigning  design  parameters  for  adaptive  optical  systems,  which  can  greatly  reduce  the 
effect  of  turbulence. 

This  reportt  describes  our  radiosonde  C,2  model*  and  tests  of  It.  Note  that  this  model 
does  not  relate  to  the  convective  boundary  layer  and  our  tests  will  be  applied  with  this  in 
mind.  Other  models  exist  for  the  boundary  layer  and  these  would  be  added,  for  example,  when 
parameters  such  as  rG,  the  coherence  length,  are  estimated,  because  ro  Is  sensitive  to  near -ground 
Cn2.  In  this  report  we  will  only  consider  parameters  sensitive  to  Cn2  above  the  boundary  layer. 

This  report  describes  how  the  AFGL  model  was  created  from  very  high  resolution  veloc¬ 
ity  profiles  that  we  obtained  in  the  stratosphere  by  means  of  rocket  laid  smoke  trails;  next  It 
will  explain  how  the  resulting  model  is  used  to  convert  radiosonde  data  Into  Cn2  profiles.  Then, 
comparisons  between  model  and  thermosonde  profiles  of  Cn2  will  be  made.  In  addition,  the 
AFGL  model  will  be  compared  to  two  other  radiosonde  models,  namely,  those  due  to  Hufnagel3 
(HUF)  and  to  VanZandt-’  (VZ).  Three  field  programs  will  be  discussed,  namely,  CLEAR- 1,  and 
Hawaii  (1984)  In  the  main  text,  and  CLEAR-3  In  Appendix  1. 

At  the  end  of  the  report  sonv-  comments  will  be  made  about  the  lessons  we  learned  In 
this  research.  Suggestions  for  future  research  will  also  be  offered.  Appendix  A  contains  a  paper 
that  describes  the  VZ  and  HUF  models  and  additional  model  comparisons,  and  Appendix  B 
contains  more  information  on  high  resolution  shear  regressions. 

Received  for  publication  26  Feb  1993 
tThe  original  draft  of  tins  manuscript  originated  25  Feb  1987. 

*  Hereafter  this  will  be  denoted  Hie  AFGL  model,  for  reasons  of  established  conventions  in  the  literature. 


2.  The  AFGL  Radiosonde  Model  for  Cn2 
(Above  Boundary  Layer) 

2.1  Basic  Concepts 

Till*  key  equation  lor  the  AFGL  model  Is: 


where: 


G=2  8A/  /; 
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and  where  T  Is  absolute  atmospheric  temperature  In  °K.  P  is  pressure  In  mb,  y  is  the  dry 
adiabatic  lapse  rate  ot  9.8  \  10  3  °K/m,  and  z  is  the  height  above  ground.  (Tatarski,  1961. 
and  Dewan,  1980).  Radiosondes  give  us  P  and  T  directly,  but  L  in  Eq.  (1)  is  “the  outer  length." 
that  is.  the  largest  scale  ot  inertial  range  turbulence.  This  is  the  unknown  that  our  model 
will  supply.  A  good  rule  of  thumb,  based  on  information  in  Pend  el  al  (1963)  and  Tennekes 
and  Lumley  (1972),  is  that  L  would  be  of  the  order  of  01  times  the  thickness  of  a  turbulent 
layer.  In  principle  this  could  be  used  in  the  future  as  an  adjustable  parameter  (that  is,  one 
could  use  values  other  than  0.1)  but  at  present  using  other  values  seems  to  be  unnecessary. 

We  now  consider  the  question  "How  can  one  estimate  L  from  radiosonde  data?"  which 
will  occupy  us  until  Section  2.3.  It  Is  generally  known  that  above  the  convective  boundary 
layer,  atmospheric  turbulence  occurs  In  thin  layers  shaped  like  pancakes  that  are  miles  in 
width  and  tens  of  meters  thick  (usually).  A  shear  type  of  instability  leads  to  the  formation  of 
the  layers,  and  the  shears  are  generally  caused  by  gra*ity  waves  (Dewan  and  Good,  1986). 7 
We  used  the  common  rule  of  thumb 


Ri  =  iV'  /  S:  <0.25 


(3) 


where  Ri  Is  the  Richardson  number,  N  is  the  buoyancy  frequency,  and  S  is  the  vector  vertic  il 
shear  of  the  horizontal  velocity  defined  as 


S  = 


dV,\2 

dz. 


+ 


(4) 


where  VN  and  VE  are  the  north  and  east  horizontal  wind  components.  In  a  paper  as  yet  un¬ 
published,  by  Dewan  and  Good,  we  found  that  data  supported  the  instability  criterion 

(S.)2  >{N)Z  (0.5)' 1  (5) 


and,  in  the  present  report  we  shall  use  Eq.  (5). 

Standard  radiosondes  report  data,  such  as  velocity,  at  Intervals  of  300  m  and  larger.  When 
such  velocities  are  used  in  Eq.  (4),  only  rarely  will  condition  (5)  hold.  Presumably,  since  one 


expects  that  the  layers  are  of  order  1/10  the  resolution  ol  the  radiosondes,  the  shears  re¬ 
sponsible  must  be  on  that  same  scale.  Van  Zandt  ct  al.  (19b l)1  pointed  this  out  and  indl 
cated  tile  resultant  need  tor  a  statistical  model  to  estimate  small  scale  wind  struct  are.  We 
must  use  a  statistical  association  between  the  large  scale  shears,  called  Sraw  here  (that  are 
measured  by  radiosondes)  and  the  average  of  L4/;i  contained  within  the  300  m  height  range, 
which  are  based  on  ttie  small  scale  but  unmeasurable  shears. 

2.2  The  Small  Scale  L-Model 

To  obtain  the  stattstieal  association  mentioned  above  we  used  our  high  resolution  (10  m) 
stratospheric  velocity  profiles  described  In  Dewan  et  al.(1984).8  In  all.  these  data  consisted  of 
55.3  km  of  velocity  profile  information.  The  choice  of  300  in  for  "radiosonde  scale".  Is  not  as 
arbitrary  as  might  appear.  It  is  based  in  part  on  subsequent  model  performance. 

Our  procedure  ts  described  with  reference  to  Figure  1 .  The  first  step  is  to  obtain  the  high 
resolution  shears  (from  the  previously  mentioned  10  m  winds)  as  a  function  of  z  by  using  Eq. 
(4).  bsing  Eq.  (5)  and  the  US  Standaru  Atmosphere  we  found  that  Sc  =  0.015  s  '  for  the  tro¬ 
posphere  and  0.03  s  1  for  the  stratosphere.  Next,  all  shear  regions  exceeding  Sc  are  presumed 
turbulent.  As  shown  In  Figure  1,  Mrbulent  layer  thicknesses,  ;f,  are  assigned  to  these  re¬ 
gions.  We  will  use: 


/,  =  1/102,  (6) 

but  the  1/10  factor  is  not  taken  into  account  until  the  model  Is  actually  applied  to  the  radio 
sonde  data. 

Next,  we  obtained  a  weighted  average  of  (./.', |4/l  from 

<^,,>S£^:/,a’,/290)  (7) 

Occasionally  a  layer  L,  will  protrude  outside  of  the  300  m  region  of  Interest.  In  such  a 
case  the  In  the  term  (i/,/2.90)  of  Eq.  (7|  Is  reduced  to  include  only  the  amount  of  if,  within 
the  300  m  region. 

Several  methods  to  obtain  SroW  are  described  in  Appendix  13;  however,  in  the  model  of 
choice,  (on  the  basis  of  performance)  we  first  smoothed  the  velocities  with  an  11  point  run¬ 
ning  average  and  then  differences  the  velocities  across  the  300  m  region  and  divided  the  re¬ 
sult  by  300  in.  This  gave  us  S.aw.  Finally  wc  plotted  <if,4/3>  against  Sraw  in  Figures  2  and  3 
on  a  log-linear  plot  and  obtained  a  straight  line  regression  In  these  coordinates.  The  piaxl- 
mum  raw  shear  used  in  the  regression  was  0.045  s  l.  The  resulting  models  are: 

log  <  X']n  >=Y  =  1 .64  +  42.0  SIilH.  (Troposphere) 

Y  =  0.506  +  50.0  Sww  (Stratosphere)  (8) 


This  Is  to  be  used  In  Eq.  (10)  below, 
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Figure  1  —  Method  for  assigning  Sriw  to  (L4/3). 
(Uses  weighted  average,  Eq  (7).) 


w((n4/')-K 


TROPOSmiSRl- 


I 

^raw  ^  Cutoff 


Figure  2  —  10T  vs  SrBW  Regression.  This  is  the  Regression  for  the 
Tropospheric  Model  [Sc  =  0.015). 
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These  regressions  were  obtained  from  data  analyzed  for  every  10  m  (that  is,  the  300  m 
region  was  shifted  by  10  m  for  each  output  pair  of  <Ll4/3>  and  Sraw).  It  is  clear  that  the  num¬ 
ber  of  Independent  points  was  1/30  of  the  total  number  of  points  (the  latter  being  of  order 
4.000).  For  this  reason  only  every  30th  point  is  exhibited  in  Figures  2  and  3.  Taking  this  into 
account,  we  estimated  (slightly  overestimated)  the  standard  deviations  for  the  regressions  as: 

OW,  “6.9 

=6x10" 


These  are  based  on  (Bevington  (1969)) 

[l/(/V-2)]£(y,-u-/;a-)2 

<v/5>.>2 

°r«2/N 

Using  a  representative  value  of  Sriiw  =  0.02  s  1  we  estimate  that  the  spread  in  the  measure¬ 
ments  introduces  an  uncertainty  of  about  a  factor  of  1 .6  in  Cn2.  Generally  such  a  spread  is 
unsatisfactory;  but  since  Cn2  can  vary  over  a  range  of  several  orders  of  magnitude,  a  factor  of 
1.6  is  considered  to  be  essentially  negligible.* 

Figures.  2  and  3  show  a  large  scatter  of  data  points  about  the  lines.  In  spite  of  this  we 
will  see  that  the  model  performs  well.  The  explanation  lies  in  the  fact  that  the  standard  de¬ 
viation  of  a  linear  regression  Is  like  that  of  a  mean  quantity,  an  additional  factor  of  the  square 
root  of  a  number  of  independent  cases  is  involved  in  the  "standard  error"  of  the  mean. 

While  Eq.  (8)  Is  the  model  of  choice,  a  slightly  different  model  was  employed  below,  namely: 


<7 


f>‘ 


slope 


or 


troposphere  Y  -  1 .57  +  40.0  Sraw 
stratosphere  Y  —  0.503  i  5 1 .2  Sra 


(9) 


but  the  difference  Is  not  Important  to  model  performance. 

An  additional  point  must  be  made  regarding  the  application  of  our  model  to  the  tropo¬ 
sphere.  As  has  been  mentioned,  all  of  the  velocities  in  our  statistical  data  base  arc  from  the 
stratosphere.  For  this  reason  we  must  assume  that  this  data  base  does  not  differ  significantly 
in  shear  statistics  from  actual  tropospheric  statistical  characteristics.* 


‘  (If  Sr;lw  =  0.04  s'1,  the  uncertainty  In  Cn2  or  10Y  (due  to  the  aconst  =  0.06  and  osiope  =  6.9)  Increases  to  a 
factor  of  ahout  2.2.  Sraw  =  0.04  s'’  is  very  rarely  seen.)) 

+  (Thls  assumption  Is  weak,  and  If  future  experiments  were  to  show  that  the  tropospheric  model  is  Inferior  to  tiie 
stratospheric  model,  it  Is  precisely  here  that  one  should  start  to  look.  As  shown  In  Dewaii  and  Good  H986|  the 
velocity  spectra  are  significantly  different  In  the  troposphere  as  compared  to  the  stratosphere.) 
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2.3  Model  Application  To  Radiosonde  Data 

Our  model  consists  essentially  of  Eq.  (1).  Equation  (2)  is  evaluated  directly  irom  raw 
radiosonde  Information.  When  we  used  high  resolution  pressure  and  temperature  sensors, 
we  preprocessed  the  data  in  a  manner  to  be  described.  Using  Eq.  (4),  the  velocity  information 
was  converted  into  shears,  which  were  then  inserted  into  Eq.  (9)  to  obtain  Y  s  log10  <i£4/3>. 
Finally,  the  C,2  (z)  profile  was  obtained  from: 

C^(z) -2.S  (0.1)4/1  Atf"  10*  (10) 

where  the  factor  (0.1)  comes  from  Eq.  (6).  Note  that  Eq.  (9)  contains  in  effect  two  independent 
pieces  of  the  model  that  are  applied  to  the  stratosphere  and  troposphere  separately.  Thus,  to 
apply  this  AFGL  model  one  must  first  ascertain  the  altitude  of  the  tropopause.  In  all  the  cases 
we  studied  this  was  unambiguous:  but  this  may  or  may  not  be  a  problem  when  the  model  is 
applied  to  a  larger  volume  of  data,  since  tropopause  height  can  be  ambiguous,  as  is  known 
from  published  data.  In  connecting  the  tropospheric  model  to  the  stratospheric  model  we  used 
linear  interpolation.  Finally,  it  should  be  mentioned  that  it  was  necessary  to  build  certain 
upper  limits  Into  the  model.  If  the  raw  shear  is  greater  than  0.04,  the  model  assigns  the  value 
0.04  to  the  shear.  This  was  necessary  because  a  single  outlier  shear  due  to  an  artifact  could 
totally  dominate  the  effects  of  the  entire  Cn2  profile. 

Similarly,  we  also  placed  an  upper  limit  on  dT/dz  in  the  tropospheric  model.  Y/e  never 
allowed  it  to  exceed  zero  and  all  gradients  which  did  were  set  equal  to  10  s.  The  stratospheric 
model  would  most  likely  require  a  positive  upper  limit  on  the  temperature  gradient  but  this 
is  left  to  future  research  because  raw  radiosonde  information  in  that  region  was  limited  when 
this  report  was  written. 


Q  -  *  -  i  m.  _j 

o.  inouei  lesis 


As  was  mentioned  in  the  introduction,  we  compared  the  model  Cn2  profiles  to  both  ex¬ 
perimental  profiles  and  to  the  profiles  generated  by  the  VanZandt4  and  Hufnagel3  models.  These 
latter  models  are  described  In  Appendix  A.  To  parameterize  the  relative  performance  of  these 
models,  we  constructed  two  “figures  of  merit”  based  upon  two  crucial  adaptive  optical  param¬ 
eters.  These  were  the  isoplanatic  angle,  0O2,*  and  the  variance  due  to  scintillation,  ox2.  These 
are  obtained  from: 

e»,w,  =  {2-9lCj'‘C.;V)U)''Vz} (11) 


•Gy  is  the  angle  over  which  ttic  wave  front  correction  (in  an  adaptive  optics  system)  Is  essentially  homogeneous.  In 
other  words,  It  Is  the  field  angle  over  which  the  optical  transfer  function  Is  constant. 
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where  k  Is  the  optical  wave  number  —2 tt/A,  and 

=  <ln{Aj A,)3  >  =  0.56A:7'f!  j Cr  ( z)(z)V(,<fe 

We  thus  use,  for  figures  of  merit, 

Z,,,sr“c;(z)U),/,dz 

*5  km 

and 


tz  max 

*L  c-'<z>«  '* 


(12) 


(13) 


(14) 


The  lower  limits  on  these  integrals  were  chosen  to  avoid  the  boundary  layer,  and  all  altitudes, 
z,  are  distances  above  ground  as  opposed  to  "above  sea  level". 


4.  Comparisons:  Clear- 1  Data 

During  the  CLEAR- 1  program  we  obtained  data  from  49  thermosonde  flights.  To  com¬ 
pare  results  from  standard  radiosondes  we  used  the  routine  meteorological  flights  from 
Paso  located  about  50  miles  from  the  thermosonde  launch  site.  Only  those  thermosonde  flights 
nearest  in  time  to  radiosonde  launches  were  used  in  our  comparisons.  Table  1  lists  the  flight 
Information.  Figures  4,  5,  and  6  show  examples  of  the  profiles.  Tables  2  and  3  list  the  figure 
of  merit  comparisons,  and  Table  4  gives  the  percentage  error  comparisons  of  Z5/3 


Table  1.  CLEAR  !  Flights 


Flight 

Number 

1V484 

Date 

Launch,  PM 

Thermosonde 

Serial  No. 

El  Paso 

Thermosonde 
Time  (Local,  S.T.) 

Radiosonde 

Local  S.T. 

Altitude 
Ranges  (km) 

1 

9/4 

0508 

1:48 

5:00 

b.3  - 

Id. 3 

2 

9/5 

0506 

1:05 

5:00 

6.0  - 

15.3 

3 

9/6 

1885 

2:05 

5:00 

5.7  - 

15.3 

4 

9/7 

1882 

2:17 

5:00 

6.1  - 

15.2 

5 

9/8 

0524 

1:43 

5:00 

5.8  - 

15.3 

6 

9/10 

1878 

5:40 

5:00 

6.1  - 

15.1 

7 

9/11 

1889 

8:03 

5:00 

8.8  - 

15.2 

8 

9/13 

2703 

2:01 

5:00 

5.9  - 

15.3 

9 

9/16 

2757 

7:51 

5:00 

5.2  - 

15.2 

10 

9/18 

2769 

1:53 

5:00 

5.3  - 

15.4 

11 

9/24 

2762 

8:08 

5:00 

5.2  - 

15.2 

(Launch) 

(TROPO  ) 
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Table  2.  Isoplanatic  Angle  Figures  of  Merit  "Z5/3” 


Radiosonde  Models 

Thermosonde 

Flight 

“Z5/3“  x  107t 

Z573“  x  107 

Number 

AFGL 

HUF 

VZ 

(Experiment) 

1 

3.93* 

1.84 

3.35 

4.50 

2 

8.27* 

1.30 

3.52 

10.9 

3 

4.63 

1.49 

3.40* 

3.81 

4 

7.83 

3.36* 

13.0 

2.67 

5 

7.12 

2.80* 

5.92 

3.76 

6 

6.62  (71)** 

1.09 

2.25* 

2.29 

7 

4.46 

0.618 

2.39* 

2.42 

8 

4.84  (6)** 

1.25 

2.42* 

2.45 

9 

4.43* 

0.665 

2.72 

5.05 

10 

7.13 

0.860 

3.83* 

4.79 

11 

7.50  (9)** 

4.93* 

5.42 

2.03 

TOTAL 

3 

3 

5 

•  =  closest  model 

t  MODEL.  1.566 

4  40.03  Sraw  (TROPO.  almost  exclusively)  0.5033  + 

51  32  Sraw  (STRAT.)  Only  two  significant 

figures  are  really  involved  here,  as  well  as  in  all  the  other  tables  below.  These  two  models  are  in  reference  to 

Eq.  (81  for  Y. 

Aii  liuiViuci^  »n 

1 1 _ _  .  _  t_  1 _ _ _  _ _  _ _ 

mcr»c  luuica  m  paioi 

itiicscs  sre  bused  on  the  utiSu 

turated  mode!  when  thl 

c  differs  from  the 

saturated  model.  Note  that  they  occasionally  are  "wild."  thus  the  need  for  saturation, 
f  This  means  that  Z5/3  Is  first  multiplied  by  107  and  the  result  Is  entered  into  the  table. 

Table  3.  Scintillation  Variance  Figures  of  Merit  “Z5/«” 

Radiosonde  Models 

Thermosonde 

Flight 

“Z5/6”  x  10 10 

“Z5/6"  x  1010 

Number 

AFGL 

HUF 

VZ 

(Experiment) 

1 

1 

i  n  -7+ 

i  .  i  • 

0.857 

1  7  1 

A  .  •  4 

1.78 

2 

4.25* 

0.595 

1.73 

4.10 

3 

2.41  (3) 

0.688 

1.84* 

1.73 

4 

3.14 

1.48* 

5.63 

1.14 

5 

2.8 

1.25* 

2.82 

1.57 

6 

2.70  (30) 

0.502 

1.05* 

1.08 

7 

1.64 

0.252 

0.982* 

0.918 

8 

2.42 

0.575 

1.20* 

0.986 

9 

2.05* 

0.348 

1.33 

2.36 

10 

3.38 

0.432 

1.85* 

2.35 

1 1 

3.70 

2.21* 

2.73 

0.773 

TOTAL 

3 

3 

5 

*  =  closest  model 

MODEL:  1.566  4  40.03  Sraw  (TROPO.  almost  exclusively)  0.503  r  51.32  Smw  (STRAT.) 
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Table  4.  Percentage  Error  Comparisons  for  Isoplanatic  Angle 
{(Model  -  Exper.)  /  Exper.]  x  100%) 


Flight 

Number  AFGL  HUF  VZ  (Zs/3  ONLY) 


1 

-13 

-59 

-25 

2 

-24 

-88 

-67 

3 

22 

60.9 

10.8 

4 

193 

25.8 

386 

5 

89 

-25.5 

57.5 

6 

189 

-52.4 

-1.75 

7 

84 

-74 

-1.24 

8 

98 

-49 

-1.22 

9 

-12.3 

-86 

-46 

to 

49 

-82 

-20 

1 1 

269 

142 

160 

Up  until  this  point  our  conclusions  (based  on  Tables  1-4)  are  that  (a)  the  VanZandt  model 
gave  the  closest  values  to  the  thermosonde  the  largest  number  of  times,  and  (b)  It  also  came 
the  closest  In  value  If  one  eliminates  the  two  worst  cases  In  all  of  them. 

Tables  2-4  were  supplemented  by  Tables  5  and  6,  which  used  an  alternative  model  as 
Indicated.  Despite  the  large  apparent  difference  in  these  regressions,  as  seen  in  Figures  2 
and  3  In  comparison  to  Figure  7  and  in  the  models,  as  Indicated  explicitly  in  the  tables,  the 
results  were  hardly  affected,  and  the  relatb'e  scores  of  the  models  did  not  change.  Figure  7 
shows  the  alternative  tropospheric  model  (the  one  used  almost  exclusively  In  the  sense  that 
the  stratospheric  model  was  not  used  significantly,  since  CLEAR- 1  radiosonde  data  were  al¬ 
most  entirely  limited  to  the  tropospheric  altitudes).  Over  the  Srav/’s  of  interest,  the  regression 
almost  overlays  the  one  in  Figure  2  and  this  explains  why  the  results  were  "robust”  to  model 
change. 
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Table  5.  Isoplanatlc  Angle  Figures  of  Merit  From  Alternative  Model  “Z5/3” 


Flight 

Number 

AP'GL 

Radiosonde  Models 
“Z5/3“  x  107 

HUF** 

VZ** 

Thermosonde 

“Z5/3"  x  107 
(Experiment)** 

1 

3.45* 

1.84 

3.35 

4.50 

2 

6.92* 

1.30 

3.52 

10.9 

3 

5.00 

1.49 

3.40* 

3.81 

4 

4.91 

3  36* 

13.0 

2.67 

5 

6.48 

2.80* 

5.92 

3.76 

7 

3.97 

0.618 

2.39* 

2.42 

8 

5.24 

1.25 

2.42* 

2.45 

9 

3.99* 

0.665 

2.72 

5.05 

10 

6.31 

0.860 

3.83* 

4.79 

11 

6.79 

4  93* 

5.42 

2.03 

TOTAL 

3 

3 

4 

*  =  closest  mode! 

**  exactly  as  in  Table  2 

MODEL:  1.566  +  29.62  Sraw.  Troposphere 

0.5084  -t  37.0!  Siaw.  Stratosphere 

Table  6.  Scintillation  Variance  Figures  of  Merit  From  Alternative  Model  “Z5/0" 

Flight 

Number 

AFGL 

Radiosonde  Models 
"Z5'6”  x  10'° 

HUF** 

VZ** 

Thermosonde 
“Z5/6”  x  10 10 
(Experiment) 

\ 

1 

1.55 

0.857 

1.71* 

1.78 

2 

3.51* 

0.595 

1.73 

4.10 

Q 

V7 

2.74 

0,688 

1.84* 

1.73 

4 

2.00 

1.48* 

5.63 

1.14 

5 

2.50 

1.25* 

2.82 

1.57 

7 

1.40 

0.252 

0.982* 

0.918 

8 

2.15 

1.575 

1.20* 

0.986 

9 

1.85* 

0.348 

1.33 

2.36 

10 

2.97* 

0.432 

1.85 

2.35 

1 1 

2.73 

2.21* 

2.75 

0.773 

TOTAL 

3 

3 

4 

*  =  closest  model 
•*  exactly  as  in  Table  3 

MODEL:  1.566  +  29.62  Sraw,  Troposphere 
0.5084  +  37.01  Sraw.  Stratosphere 


NOTE:  In  Tables  5  and  6,  flight  number  6  was  omitted. 
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KILOMETERS 


KILOMETERS 


AFGL  ESTIMATE  Cn**2  1IUFNAGEL  ESTIMATE  Cn**2 


VAN  ZANDT  ESTIMATE  CV"2 


Figure  5  —  #4  in  Tables  1,  2,  and  3.  (Model  and  Data  Comparisons  in  CLEAR-1). 


KILOMETERS 


nr11  it) nr1"  10  '*  nr17  io 10  "  10  14  m1' 
'  VAN  ZANDT  ESTIMATE  Cn**2 


Figure  6  —  #10  in  Tables  1,  2,  and  3.  (Model  and  Data  Comparisons  in  CLEAR- 1) 


5.  Comparisons,  Hawaiian  Program  1984  Data 


In  1984  we  conducted  a  Held  program  In  Hawaii  during  which  a  number  of  Instrument 
packages  containing  1’  and  T  sensors  were  launched  and  lor  which,  a  tracking  system  gave 
velocities.  These  data  differed  in  several  ways  from  standard  radiosonde  data  but  they  never¬ 
theless  provided  useful  Information.  They  gave,  for  example,  an  opportunity  to  perform  strato 
spheric  comparisons  (an  opportunity  not  given  In  CLEAR- 1).  Secondly,  the  P  and  T  measure¬ 
ments  were  provided  at  a  relatively  high  resolution  (around  20  m)  in  comparison  to  radio¬ 
sondes  (300  m  and  larger).  In  this  way  we  could  better  ascertain  correlations  between  model 
and  <ri  situ  measurements  of  Cn2. 

Again  we  employed  the  model  given  In  Eq.  (9);  but  this  time  the  P  and  T  data  were  pre- 
processed.  Figure  8  shows  an  example  of  raw  P,  T.  and  velocity  data.  The  velocity  resolution 
(unlike  P  and  T)  was  very  much  like  that  found  In  standard  radiosondes.  For  this  reason  the 
velocity  data  were  not  pre-processed. 

To  smooth  P  and  T  (and  report  out  M2  at  equally  spaced  altitudes  of  300  in)  we  fitted  an 
18  point  parabola  around  the  altitude  of  interest  (that  Is,  at  each  extreme  of  the  300  m  Inter¬ 
val  for  T).  The  results  were  reported  at  the  midpoints  of  these  Intervals  for  AT/ A/  and  then  T 
and  P  were  calculated  (in  the  same  manner)  at  those  same  midpoints.  Then  the  shears  were 
calculated  from  the  raw  velocities  and  the  mode)  |Eq.  (10))  was  applied  as  was  described  In 
Section  3.  Figures  9  and  10  give  the  two  cases  that  extended  into  the  stratosphere  (the  oth¬ 
ers  did  not).  The  correlations  were  encouraging.  Tables  7  and  8  are  based  on  the  figure  of 
merit.  In  this  case  the  AFGL  model  performed  better  than  the  other  models,  but  raw  radio 
sonde  data  represents  the  type  of  data  tt  must  use  In  our  application.  Nevertheless  these  results 
argue  for  more  research  to  be  done,  in  view  of  the  success  of  the  AFGL  model. 

One  detail  regarding  Table  7  should  be  mentioned.  When  It  was  constructed  (1985)  we 
did  not  place  the  lower  limit  of  the  integral  at  5  km  but  at  the  altitude  where  data  commenced. 
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Minimum  Z  =  2.969768  Maximum  Z  =  26.595832 

Minimum  Temp  —  192-393163  Maximum  Temp  =  2!'3. 643674 


ALT1TUC 


Minimum  Z  =  2.969768  Maximum  Z  =  7.6.595832 

Minimum  V  N-S  =  9.677270  Maximum  V  N-S  =  4.823946 


Minimum  Z  =  2.969768  Maximum  Z  =  26.595832 

Minimum  Press  =  20.768661  Maximum  Press  =  706.637421 
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Figure  9-A  —  Models  and  Data  Comparisons  KAW  01  (1984). 


THERMOSONDE 


TTmiii 


Tinnir 


TTT 


"TTT 


tti  n  in  r 


i  i  i  rnnr 


1  I  I  llll 


10 


-20 


10 


-IS' 


10 


-18 


10 


-17 


10 


-If) 


10 


-IS 


Figure  9-B  —  Models  and  Data  Comparisons  HAW  01  (1984). 
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Figure  10-A  —  Models  end  Data  Comparisons  HAW  49  (1984). 
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Table  7.  Isoplanatic  Angle  Figures  of  Merit  for  the  1 984  Hawaiian  Program  Data 

[(ZB/3)  3/5  x  10  31 


RUN 

AFGL 

Thcrmosondc 

HUF** 

VZ** 

(Experiment) 

1 

6.72* 

8.91 

8.89 

7.77 

22 

1  1.6* 

37.3 

16.0 

1  1.8 

37 

3.33* 

5.66 

1.88 

3.53 

38 

6.80* 

10.2 

4.48 

6.20 

45 

7.86 

14.5* 

4.08 

12.4 

49 

5.69* 

6.42 

7.51 

4.00 

TABLE  8. 

Scintillation  Variance  Figures  of  Merit  for  the  1984  Hawaiian  Program 

Data  {(ZB/a)  x  1010] 

T*»l 

r\u  l\ 

AFGL 

Thermosonde 

HUF** 

vz** 

(Experiment) 

i 

2.41 

1.47* 

1.62 

1.56 

22 

1.56* 

0.258 

0.877 

1.56 

37 

6.15* 

3.11 

12.8 

6.19 

38 

2.43* 

1.28 

4.36 

2.28 

45 

2.62* 

0.861 

6.95 

1.33 

49 

3.24* 

2.47 

2.83 

4.93 
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6.  Comments,  and  Suggestions  for  Future  Work 


If  one  assumes,  as  we  have,  that  the  only  useful  application  of  radiosonde  type  models 
is  in  connection  with  standardized  radiosonde  data,  then  it  follows  that  future  research  should 
use  such  standardized  data  exclusively.  It  may  be  objected  that  high  resolution  data  obtained 
from  the  thermosonde  balloon  package  are  preferable  because  of  the  advantage  of  colloca¬ 
tion;  and,  furthermore,  one  could  smooth  this  data  to  a  resolution  of  300  m  to  simulate  ra¬ 
diosonde  data.  Our  experience  directly  contradicts  this  and  we  therefore  strongly  recommend 
against  further  work  with  non-standard  radiosonde  data.  Further  arguments  can  be  supplied; 
but,  a  model  should  be  designed  for  the  data  intended. 

One  possible  extension  of  the  AFGL  model  that  we  did  not  test  was  a  model  based  upon 
a  variable  Sc  value.  This  may  be  worth  pursuing  in  the  future.  At  this  point  we  would  like  to 
mention  some  of  the  attempts  that  did  not  seem  to  be  beneficial  in  this  projeci.  This  will  help 
others  who  might  want  to  continue  the  work  to  avoid  the  same  pitfalls,  At  one  point  we  de¬ 
signed  a  model  with  a  variable  resolution,  (for  example,  300  m,  600  m,  1  km,  etc.)  and  let  the 
actual  resolution  pick  the  Y  regression  to  be  used.  These  variable  models  consistently  per¬ 
formed  less  effectively  than  our  fixed  300  m  model.  Perhaps  our  high  resolution  velocity  data 
base  was  too  small  and  perhaps  a  future  effort  along  these  lines  should  be  based  on  such 
scaling  laws  as  are  found  in  the  work  of  Essenwanger. 10  Another  unfruitful  avenue  was  the 
attempt  to  improve  the  Y-regressicn  by  using  a  parabola  instead  of  a  straight  line.  This  led  to 
bad  results  and  should  therefore  be  avoided. 

The  last  blind  alley  we  will  mention  is  our  attempts  to  use  linear  prediction  techniques 
in  place  of  the  statistical  shear  model  for  Y.  We  had  hoped  to  produce  a  linear  predictor  to 
Interpolate  shears  in  the  low  resolution  radiosonde  velocity  measurements.  This  linear  pre¬ 
dictor  was  based  upon  the  high  resolution  data.  It  failed  for  reasons  unknown  to  us  at  the 
time.  We  now  know  that  the  explanation  lies  in  the  physical  basis  of  linear  predictors.  The 
latter  is  described  In  Uewan11  11985)  and  in  essence  It  ts  that  a  linear  predictor  represents  a 
dynamical  system  (or  a  filter)  stimulated  by  random  noise,  and  the  predictions  consist  of  this 
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nals  after  the  noise  Input  Is  terminated  but  eventually  they  die  out.  Unfortunately,  the  linear 
predictor  for  velocity  fluctuations  ts  highly  damped,  so  highly  damped  that  the  fluctuations 
die  out  almost  Instantly.  For  this  reason  they  were  useless  for  interpolating  high  resolution 
(10  m)  shears  between  300  m  velocity  points. 

One  final  point  should  be  made.  It  Is  that  we  still  do  not  possess  a  reasonably  large 
number  of  stratospheric  and  tropospheric  comparisons  between  the  models  based  on  radio¬ 
sondes  under  very  different  conditions.  For  example,  under  low  wind  conditions  the  Hufnagel 
model  performs  very  badly.  Clearly  the  next  step  would'be  to  remedy  this  situation. 

It  should  be  noted  that  we  learned  during  the  course  of  this  research  that  solar  heating 
of  the  thermosonde  probes  during  the  daytime  could  cause  their  measurements  to  be  in  er¬ 
ror.  Only  night  measurements  can  be  trusted  at  this  time.  (See  Reterence  12.) 


28 


7.  Conclusions 


The  AFGL  model  Is  much  simpler  than  the  VanZandt  model.  With  the  data  available  to 
us,  it  was  not  possible  to  establish  which  of  the  two  superior  models  (compared  to  the  Hufnagcl 
model)  was  actually  the  best  on  the  basis  of  performance.  In  CLEAR- 1,  the  NOAA  (VanZandt 
and  Warnock)  model  was  rated  best;  when  the  Hawaiian  campaign  data  were  used,  the  AFGL 
model  performed  best.  If  future  tests  show  that  the  models  remain  more  or  less  comparable 
in  performance,  then  the  decision  of  which  one  to  use  could  be  made  on  the  basis  of  slmplic 
ity;  the  AFGL  model  would  then  be  chosen. 
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Appendix  A 


This  appendix  contains  the  entire  report  presented  in  1986  at  a  conference  (compare  with 
Good  et  al.  1988). A1  it  Includes  a  description  of  the  Van  Zandt  and  Hufnagel  models.  It  also 
provides  additional  model  comparisons.  Perhaps  one  more  asset  is  that  it  will  provide  an  al¬ 
ternative  description  to  compare  with  the  main  text. 


MODELING  CN2  PROFILES  FROM  RADIOSONDE  DATA 

E.  M.  Dewan,  R.  E.  Good,  R.  Beland,  J.  Brown 
Hanscom  AFB,  Bedford,  MA  01731-5000 


ABSTRACT 

The  structure  constant  for  the  fluctuations  ol  the  optical  index  of  refraction.  CN2,  is  the 
key  to  the  design  of  adaptive  optical  systems  so  that  the  effects  of  turbulence  on  laser  beam 
propagation  can  be  minimized.  This  paper  compares  three  models  which  cuuvert  standard 
radiosonde  data  into  CN2  profiles.  Models  arc  compared  to  directly  measured  values  of  CN2 
obtained  by  balloon  borne  thermosondcs. 


AX.  INTRODUCTION 

The  purpose  of  this  section  is  to  compare  the  performance  of  the  three  existing  models 
for  converting  radiosonde  data  into  CN2  vertical  profiles.  The  outputs  of  these  three  models 
will  be  compared  to  actual  in  situ  experimental  thermosonde  measurements.  In  addition  to 
direct  profile  comparison  a  parametric  comparison  will  be  made  based  on  isoplanatic  angle 
calculations. 

During  CLEAR  3  we  obtained  three  thermosonde  and  four  radiosonde  profiles  on  two 
simultaneous  dates  and  they  were  all  launched  from  WSMR.  They  were  not  launched  on  the 
same  hours  nor  from  the  same  locations  however.  Figure  Al  shows  a  comparison  of  a  radio 
sonde  temperature  profile  taken  at  one  location  and  a  thermosonde  flight  temperature  profile 
taken  on  the  same  day  but  from  a  different  location.  As  can  be  seen,  these  profiles  are  quite 
similar,  and  the  main  differences  could  perhaps  be  attributed  to  differences  in  resolution. 

In  order  to  eliminate  effects  of  boundary  layer  turbulence  (which  are  outside  of  these 
model  capabilities)  we  Intentionally  ignored  all  information  below  5  km.  This  has  the  added 
benefit  of  eliminating  the  need  to  include  humidity  effects  which,  strictly  speaking,  do  not 
extend  beyond  about  7  km. 
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Our  conclusion  will  be  that  the  results  suggest  that  the  AFGL  model  may  be  superior  In 
pertormance  and  simplicity  to  the  Van  Zandt  model.  More  data  would  be  needed  however  to 
confirm  that  this  is  true  in  general.  Also,  as  was  expected,  the  Hufnagel  model  falls  under 
the  experimental  conditions  of  CLEAR  3. 


A2.  CN2  Profile  Models 

Historically,  the  oldest  and  simplest  of  these  models  is  the  one  due  tc  Hufnagel  (Refer¬ 
ence  Al).  It  Is  given  by 

C-l  ( z )  =  [8.2  x  1CT56  V2  Z10  je‘u/10'K'1  +  (2.7  x  1CT1g  e  ~ lz/15000)  j  ( i) 

where  Z  is  the  altitude  In  meters.  V2  is  tfe  spatial  average  of  the  wind  speed  squared  over 
the  altitude  regions  on  5km  to  20km,  and  CN2  is  given  in  units  if  nr2/3.  The  simplicity  of  this 
model  resides  In  the  fact  that  only  V2,  a  single  parameter.  Is  required  of  the  radiosonde  data. 

The  second  model  Is  due  to  Van  Zandt  et  al,  (References  A3.  A4).  It  is  the  most  complex 
of  three  under  discussion.  In  order  to  assign  a  value  of  CN2  to  a  specific  height  region  one 
needs  a  number  of  parameters  from  the  radiosonde  data.  These  arc  temperature,  T,  in  ‘K, 
pressure  P  in  mb,  temperature  gradient.  dT/dz.  specific  humidity,  q.  and  its  gradient  dq/dz, 
and  the  vector  wind  shear  from  the  velocity  measurements,  S.  Since  this  model  is  based  In 
part  upon  the  Richardson  number  (Ri)  criterion,  which  is  given  by  Ri  =  N: /S2  <  1/4.  where  N 
is  the  buoyancy  frequency,  the  two  key  parameters  are  given  by  N2  and  S2.  That  is  to  say  N2 
and  S2  must  be  fed  Into  the  n.odel  on  the  basis  of  the  data.  As  has  been  mentioned  already, 
we  will  omit  the  use  of  q  and  dq/dz  In  this  report  (q  Is  specific  humidity). 

A  second  and  crucial  piece  cf  this  model  is  given  by  the  relation 

C'N-'lX>  M'  l'  (2) 


(See  Tatarskl,  Reference  A5)  where  L  j  the  "outer  length".  In  the  present  context  this  is  taken 
to  mean  a  representative  scale  associated  with  turbulent  layer  thickness.  The  parameter  M. 
which  is  the  factor  giving  the  gradie''*  of  the  index  of  refraction,  and  hence  the  bridge  be¬ 
tween  L  and  CN2,  is  given  by 

..  -79x10  -*PN2 


where  P  is  in  mb  and  T  in  'K,  N  is  s  1  and  g  is  the  acceleration  of  gravity.  The  definition  of  N 
is  given  by 


A,:  _  jKf + r) 


(4) 
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where  y  ts  the  dry  adiabatic  lapse  rate  of  9.8  x  10  3  fK/ni).  In  the  above  equations  suitably 
averaged  quantities  are  assumed. 

In  the  development  of  this  model.  Van  Zandt  et  al.  assume  that  shears  and  temperature 
gradients  on  a  very  small  scale  (relative  to  the  scale  of  radiosonde  data)  are  the  cause  of  the 
turbulent  layers.  This  smaller  scale  they  refer  to  as  the  microscale.  They  make  a  number  of 
assumptions  regarding  the  probability  distribution  of  microscale  wlndshears  and  temperature 
gradients.  For  example  they  assume  a  Rtce-Nakagami  distribution  for  the  absolute  values  of 
vector  wind  shears  and  a  Gaussian  distribution  for  the  temperature  gradients  affecting  N2. 
They  also  assume  that  the  Joint  probability  distribution  for  L,  S2,  N2,  and  dq/dz  can  be  fac¬ 
tored.  They  also  made  use  of  the  statistical  findings  of  Essenwanger,  Rosenberg  and  Dev/an, 
Cadet,  and  Adel  fang  (see  References  A3  and  A4  for  the  references).  Finally  their  model  is  predi¬ 
cated  on  the  assumption  that  the  turbulence  under  consideration  Is  caused  by  gravity  waves. 
Reference  3  should  be  consulted  for  all  the  details  as  well  as  for  a  very  well  documented  com¬ 
puter  program  for  this  model.  As  can  be  seen.  It  Is  apparently  somewhat  complex.  The  key 
equation  for  the  Van  Zandt  model  Is 

Cl  (z)  =  j dt\  J dSdN  P'  ( L,S,  N ) C;  ( L  N)  5(a) 

() 

The  third  model  we  shall  call  the  “AFGL  model."  It  is  the  most  recent  one  to  date  and  it 
Is  simpler  than  the  above  case  despite  the  fact-  that  the  physics  behind  it  is  essentially  the 
same  In  the  sense  that  both  Rl  and  Eq.  (2)  enter  into  its  construction.  Its  Implementation 
relies  on  the  simple  relation 


C2n  -  2.8  [f..'(z)]2  (0.1)4/3  10vu)  (5b) 

As  before,  P  and  T  from  the  radiosonde  determine  M2  from  (3)  and  (4).  and  Eq  (2)  leads  to  (5). 
The  factors  i  10Y(0. 1)4/3|  determine  L4/3  in  a  statistical  manner  to  be  described.  The  factor  (0.1) 
arises  from  a  rule  of  thumb  that  the  outer  length,  L  is  of  the  order  of  one  tenth  the  thickness 
of  a  turbulent  layer.  In  the  Implementation,  Y  is  obtained  from 

Y  =  C,  +  C2  5  (6) 

where  S  is  the  “raw  shear"  (not  the  microscale  shear)  directly  measured  from  the  radiosonde 
data.  The  term  Y  arises  from  the  definition  log  i£4/3  a  Y,  and  10Y  determines  the  most  prob¬ 
able  median  value  of  L4/3  in  the  altitude  region  under  consideration.  Again,  following  the  lead 
of  the  Van  Zandt  approach  this  model  assumes  that  microshears  cause  turbulence;  however, 
the  value  of  N2  used  in  the  Richardson  number  criterion  is  no  longer  data  dependent  In  the 
AFGL  model.  Instead,  the  value  of  N2  is  based  on  a  model  atmosphere  during  the  creation  of 
the  model.  In  the  application  of  the  model  the  temperature  data  Is  used  to  calculate  M  in  Eq. 
(5)  as  already  mentioned.  Note  that  this  Is  a  departure  from  the  Van  Zandt  model  in  that 
microscale  temperature  gradients  are  not  in  the  theoretical  picture.  On  the  other  hand,  the 
AFGL  model  does  Incorporate  radiosonde  measurements  of  temperature  in  the  M2  factor. 


Since  the  value  of  N2  used  In  the  Richardson  number  criterion  is  not  data  dependent, 
the  AKGL  model  Is  actually  composed  of  two  models,  one  for  the  stralosphere  and  one  for  the 
troposphere.  During  application,  the  results  from  these  two  models  are  simply  Joined  by  a 
straight  line  across  the  tropopause.  These  two  models  consist  of: 

Y  =  1 .566  +  29.62  S  TROPOSPHERE  (7) 

y  =  0.508  + 37.01  S  STRATOSPHERE 


Unlike  the  Van  Zandt  model,  the  AFGL  model  has  not  been  explained  in  detail  anywhere 
in  print  and  hence  It  may  be  useful  to  describe  how  we  arrived  at  Eq.  (7).  We  were  fortunate 
to  be  In  possession  of  a  unique  set  of  very  high  resolution  wind  profiles  (10  m).  These  we 
believe  represent  actual  measurements  of  the  microscale  motions  and  hence  we,  in  a  sense, 
had  direct  access  to  the  data  relevant  to  the  microscale  turbulent  layers.  On  the  basis  of  the 
previously  mentioned  standard  atmospheric  value  of  N2  we  arrived  at  0.0155  s  1  and  0.030  s 
l*  for  the  tropospheric  and  stratospheric  values  of  the  critical  mlcroshear  (that  Is,  the  shear 
necessary  to  cause  turbulence  on  the  basis  of  the  Richardson  criterion).  In  this  way  we  could 
assign  layer  thicknesses  to  the  microscale  data.  In  other  words,  a  thickness  was  assigned 
to  each  region  where  the  microshears  exceeded  the  critical  value. 

Since  the  resolution  of  the  radiosondes  Is  of  the  order  of  300m,  we  related  the  ‘raw 
shears,”  that  is  the  shears  across  300m  of  the  microscale  winds,  to  the  value  of  log  i£4/3  where 
the  ovei  bai  icpicscuts  the  weighted  average.  Thus  we  used 


T  ‘  (300m) 


(8) 


for  each  segment.  On  the  basis  of  linear  regression  we  obtained  Eq.  (7). 

In  summary,  the  AFGL  model  is  applied  to  radiosonde  data  by  inserting  the  measured 
values  of  T.  P,  and  V  (the  vector  velocity)  and/or  the  appropriate  derivative  etc.  Into  (7)  and 
(5b)  'using  (4)  and  (3).  The  results  from  (5b)  are  then  plotted  as  a  function  of  height.  This 
model  seems  to  be  simpler  In  concept  and  application  than  the  Van  Zandt  model.  We  now 
iurn  to  the  cumpaiisun  of  results. 


*  USA  Standard  ATM  actually  Rives  data  implytnR  .025  s-'  for  the  troposphere  and  0.045  s'*  for  the  stratosphere  til 
previous  work.  The  above  values  were  chosen  because  they  Rave  better  performance  In  previous  comparisons. 
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A3.  Model  Profile  Comparison 


Figure  A2  shows  a  representative  comparison  between  the  model  CN2  proi.  and  the  in 
situ  measurements.  The  most  obvious  feature  of  this  comparison  is  the  fact  that  the  Hufnagel 
model  profile  is  far  too  low  in  CN2  values  as  compared  to  the  data.  In  this  respect  it  has  failed 
to  give  adequate  results.  There  is  reason  not  to  be  surprised  by  this  situation;  namely,  it  is 
well  known  that  the  Hufnagel  model  will  do  this  under  conditions  where  the  winds  are  rela¬ 
tively  weak  such  as  was  the  case  for  CLEAR  3. 

Another  salient  feature  is  the  fact  that  both  the  Van  Zandt  and  AFGL  models  give  lower 
CN2  values  than  the  data  in  the  statosphere.  One  possible  interpretation  is  that  these  models 
both  perform  less  adequately  m  the  stratosphere.  The  same  holds  for  the  Hufnagel  model. 
Another  Interpretation  (in  the  context  of  possible  solar  heating  effects  on  the  thermosonde 
probe)*  is  that  in  actual  fact  the  data  are  overly  large  due  to  the  heating,  in  other  words  it  is 
possible  that  the  data  and  not  the  models  are  note  to  blame.  Parenthetically  it  should  be 
menuoned  that  Van  Zandt's  (et  al.)  model  used  here  is  based  on  Reference  A3  not  Reference  A4. 

Our  parametric  characterization  is  based  on  the  well  known  and  important  isoplanatic 
angle  0O  It  is  related  to  the  CN2  profile  through  the  relation 


where 


Q  (  r/tr1\  -r 

Vti  v - / 


0  Q5  (  —  )  1 

\  X  )  * 


a/ s 


1=  Jc£(z)(z-zfl)5/3dz 


5xlO:'m 


O) 


(10) 


and  where  meters  are  used  throughout,  zg  is  the  altitude  of  the  ground  station  (zero  here  in 
the  comparison  role),  and  X  is  the  wavelength  of  the  optical  signal.  Note  the  lower  limit  on 
the  Integral.  In  practice  it  would  start  at  z^  but  our  purpose  is  not  to  calculate  the  true  0O  but 
rather  to  compare  models  and  ignore  the  boundary  layer  as  has  been  already  explained. 

Table  Ai  lists  flight  dates  and  times  as  well  as  values  of  z  max  for  overlap.  The  compari¬ 
son  parameter  I  shows  that  the  AFGL  model  gave  the  closest  agreement  with  the  data.  Table 
A2  gives  a  similar  comparison  in  term  of  and  Table  A3  gives  the  '‘percentage"  comparisons. 
If  the  thermosonde  data  are  assumed  to  be  the  final  arbiter  of  this  comparison,  then  the  AFGL 
model  was  the  best  one.  On  the  other  hand,  unanswered  questions  remain  as  have  been  pointed 
out  in  the  text.  Furthermore  we  need  more  comparisons  under  varied  conditions  before  the 
issue  can  be  settled.  Our  modest  conclusion  is  that  If  the  AFGL  model  were  merely  to  per¬ 
forin  as  well  as  the  Van  Zandt  model,  it  would  still  be,  in  a  sense,  superior  because  it  is  a 
much  simpler  model. 


•This  possibility  is  presently  under  very  active  investigation.  NOTE  added:  it  turns  out.  as  mentioned  in  (he  main  text, 
that  daytime  ihermosonde  measurements  are  unreliable.  In  fact,  they  are  now  totally  avoided  In  current  research. 
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Table  A-l.  Values  of  I  x  107  for  Various  Radiosonde  Models  and  for 
Thermosonde  Measurements  in  the  Altitude  Range  of  5  km  to  z  max 


Radiosonde  Models* 

Thermosonde* 

7 

■■■ 

x  107  ' 

Launch 

L 

Max  Alt  In 

1  x  10' 

Time 

Ser. 

Common 

Date 

AFGL 

Hufnagel 

Van  Zandt 

EXP. 

MST 

No. 

(km) 

1  Aug  1 

8.17  am 

11.5 

0.994 

5.64 

4.93 

1:44  prn 

6474 

24.98 

2  Aug  2 

0.03  am 

11.8 

1.18 

4.00 

10.1 

1 1 :08  am 

6483 

25.16 

3  Aug  2 

0.03  am 

11.8 

1.18 

3.98 

12.8 

3:19 

5306 

24.00 

4  Aug  2 

6.00  am 

6.54 

1.12 

4.71 

10.1 

1 1:08  am 

6483 

25.16 

5  Aug  2 

6.00  am 

6.50 

1.12 

4.69 

12.8 

3:19  pm 

5306 

24.00 

6  Aug  2 

9.5  am 

5.54 

0.956 

3.61 

5.50 

11:08  am 

6483 

16.9 

7  Aug  2 

9.5  am 

5.54 

0.956 

3.61 

8.27 

3:19  pm 

5306 

16.9 

•  4  Radiosonde  profiles  —  3  models,  compared  to  3  Thermosonde  profiles. 


Table  A-2  —  Values  of  0O  for  Various  Radiosonde  Models  and 
Thermosonde  Measurements 


Model  G 

0  (p-rad) 

EXP.  60 

(p-rad) 

AFGL 

Hufnagel 

Van  Zandt 

Thermosonde 

Ser.  No. 

4.58 

19.9 

7.02 

7.61 

6474 

1 

4.51 

17.9 

8.62 

4.95 

6483 

2 

4.51 

17.9 

8.65 

4.29 

5306 

3 

6.42 

18.5 

7.82 

4.95 

6483 

4 

6.44 

18.5 

7.84 

4.29 

5306 

5 

7.09 

20.4 

9.17 

7.12 

6483 

6 

7.09 

20.4 

9.17 

5.58 

5306 

7 

Dates,  times,  and  altitude  ranges  are  exactly  as  in  Table  1  In  sequence. 
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Table  A-3  —  Comparison  Between  AFGL  and  Van  Zandt  Models 


No. 

AFGL  A  % 

VAN  ZANDT  A  % 

1 

-39.8 

-7.75 

2 

-8.89 

58.8 

3 

5.13 

-102.0 

4 

29.7 

58.0 

5 

50.1 

82.8 

6 

-0.42 

28.8 

7 

27.06 

64.3 

AVE  ABS.  VAL 

23  t  19 

57  ±  30 

A  %  = 

(MODEL  0O  - 

EXPER.0J  x  100 

EXPER.0o 
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Figure  A2  —  CN2  Proflies:  Hafuagei,  Van  Zandt,  and  AFGL  models 
compared  to  thermosonde  data. 


Appendix  B 

Shear  Regressions 


We  calculated  several  regressions  for  Y  In  our  shear  model  which  will  now  be  described. 
Our  purpose  Is  to  describe  alternate  models  (some  used  In  the  text,  and  some  not  used)  In 
the  hope  that  the  results  might  prove  useful  for  future  work. 

In  the  text  we  described,  in  connection  with  Eq.  (8),  the  technique  of  using  an  1 1  point 
smoothing  to  calculate  the  Sraw  values  across  a  300  m  altitude  range  containing  an  equal  range 
(that  is,  300  m)  of  weight  averaged  length  scales  which  gave  the  values  of  Y  to  plot  against 
Sraw  This  model  was  obtained  for  both  stratosphere  and  troposphere.  Other  models  Involved 
smoothing  over  29  points  as  well  as  no  smoothing.  All  are  listed  In  Table  B 1 . 


Table  Bl.  Alternative  Models 


A.  Troposphere 

Constant 

Slope 

1.  no  smoothing 

1.65 

39.6 

2.  1 1  pt  smooth 

1.64 

42.0 

3.  29  pt  smooth 

1.65 

47.6 

B.  Stratosphere 

1.  no  smootiiing 

0.521 

47.2 

2.  1 1  pt  smooth 

0.506 

50.0 

3.  29  pt  smooth 

0.525 

56.2 

Note  that  Increasing  the  smoothing  has  little  impact  upon  the  constant  but  that  it  In¬ 
creases  the  slope.  This  makes  sense  because,  as  smoothing  Is  increased,  it  will  reduce  Sraw 
for  a  given  value  of  associated  <if4/a>,  the  weight-averaged  if4''3. 

In  the  Initial  stages  of  our  work,  we  employed  a  model  (used  also  In  Appendix  A  and  in 
the  CLEAR- 1  study)  which  consisted  of: 

Tropo  Y  =  1 .57  +  29.6  Sruw 
Strat  Y  =  0.508  +  37.0  SraM, 

as  mentioned  In  the  texts  (see  Tables  5  and  6).  This  model  was  obtained  In  a  calculation  that 
Involved  a  150m  region  enclosed  within  the  300  m  region  and  centrally  located.  It  also  in¬ 
volved  a  weight  averaged  value  of  if,  that  consisted  of: 

{(^r)=V(x;'’(2,/290)) 

[compare  Eq.  (7)  of  the  main  text)  which  made  no  allowance  for  if,  sections  extending  outside 
of  the  volume  of  Interest.  While  we  consider  this  procedure  to  be  dubious  at  present,  we  found 
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that  In  our  CLEAR-2  work  (omitted  in  this  report  but  contained  in  part  in  the  CLEAR-2  re¬ 
port)  the  model  seemed  to  work  well.  We  omitted  the  details  here  because  we  were  not  sure 
that  the  CLEAR-2  velocities  were  reliable  in  unfiltered  form  as  we  used  them.  In  any  case, 
the  text  shows  that  performance  was  very  close  to  the  1 1  point  smooth  case  given  in  Table  Bl. 


42 


